Electronic structure of strongly correlated d-wave superconductors 
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We study the electronic structure of a strongly correlated d-wave superconducting state. Combin- 
ing a renormalized mean field theory with direct calculation of matrix elements, we obtain explicit 
analytical results for the nodal Fermi velocity, vf , the Fermi wave vector, fej?, and the momentum 
distribution, rih, as a function of hole doping in a Gutzwiller projected d-wave superconductor. We 
calculate the energy dispersion, Ek, and spectral weight of the Gutzwiller-Bogoliubov quasiparticles, 
and find that the spectral weight associated with the quasiparticle excitation at the antinodal point 
shows a non monotonic behavior as a function of doping. Results are compared to angle resolved 
photoemission spectroscopy (ARPES) of the high temperature superconductors. 
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Recent progress in angle resolved photoemission spec- 
troscopy (ARPES) of the high temperature supercon- 
ductors has led to considerable interest in the electronic 
structure of a strongly correlated d-wave superconduct- 
ing state P, 0- Experiments show a variety of inter- 
esting phenomena and it is generally agreed that strong 
electronic correlations play a dominant role in explaining 
some of the universal spectral features. Hence it would 
be desirable to obtain, e.g., for comparison with exper- 
imental results, explicit analytical results from a simple 
model of a strongly correlated d-wave superconductor. 

Thus motivated, we consider in this paper, the elec- 
tronic structure of a Gutzwiller projected supercon- 
ductor, defined by the ground state wave function, 
exp(i£)|\D f ) = cxp(iS)PG\^o) ■ The projection operator, 
Pq = Hi(l — n i] n ii)^ acting on a BCS wave function, 
|>?o)) eliminates states with double occupancies in 
The operator, exp(iS), allows for a systematic calcula- 
tion of corrections to the fully projected state, |\&) 
Gutzwiller projected states, |\&), were initially proposed 
as variational states to describe superconductivity in the 
proximity of a Mott insulating phase [H, El 0- They 
have been used recently for numerical investigations of 
the electronic structure of the high temperature super- 
conductors. Paramekanti et al., used Variational Monte 
Carlo (VMC) to study some spectral properties of a 
Gutzwiller projected d-wave superconductor [2|. Yunoki 
et al., extended the VMC technique to the direct calcu- 
lation of excited states in Jastrow-Gutzwiller wave func- 
tions 8]. The VMC technique also allows to study the co- 
existence of superconductivity with antiferromagnetism 
0, and the quasiparticle current renormalization |lfj| . 

In this paper, we will follow an alternate route and 
use a combination of renormalized mean field theory 
(RMFT) and direct calculation of matrix elements 
| 111 ll| to examine the electronic structure of a Gutzwiller 
projected d-wave superconductor. Though Gutzwiller 
projection is only an approximate method to treat the 



effects of strong correlations, the advantage lies in its 
directness and clarity. Explicit analytical expressions de- 
rived in this paper can be used to evaluate directly the 
successes as well the limitations of this approach. Us- 
ing RMFT, we determine the energy dispersion, Ek, of 
the Gutzwiller-Bogoliubov quasiparticles, and the nodal 
Fermi velocity, vf- We show that vp stays finite as the 
(Mott) insulating phase is approached. We calculate the 
spectral weight associated with the d-wave quasiparticles, 
and the momentum distribution, rife, by a direct evalua- 
tion of the relevant matrix elements 11]. We find that 
the quasiparticle weight associated with the antinodal ex- 
citation exhibits a non monotonic behavior as a function 
of doping. 

We consider the one band Hubbard model, 



H = - E *w) (4 



U ^2 n it n il 



where the hopping integrals, tuj\, connect sites i and j. 
We will restrict our attention to nearest (t), and next 
nearest (<' = —t/4) neighbor hopping. We choose an 
onsite repulsion U = 12 1; i.e., we work in the strong 
coupling regime U t.t' |l2j|. In this limit, a uni- 
tary transformation, exp(iS), yields an effective Hamil- 
tonian, H c s — exp(— iS)H exp(iS), which is block diago- 
nal and does not mix states between the lower and upper 
Hubbard bands HEl. To 0{t 2 /U), it corresponds to 
the t — J Hamiltonian with the correlated hopping term 
(three-site term), for which the Gutzwiller wave func- 
tion PqY&q) has been established as an excellent vari- 
ational ground state [(|. Retransformation of the trial 
wave function, exp(iS)Pc\^o) i then provides a system- 
atic way to study the Hubbard model using Gutzwiller 
projected states. 

Two steps are necessary to obtain explicit analytic 
expressions for the low energy properties in the strong 
coupling regime. The first is the Gutzwiller renormal- 
ization procedure, where the effects of the projection 
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Pa are taken into account by appropriate renormaliza- 
tion factors, following Hilbert space counting arguments: 

(^o\p G H cfi p G \^o) « <*oii?cfri*o) 0111113 

The next step is the realization that H c s allows for two 
types of molecular-fields: the hopping amplitude, £ r = 
So- ( c la c i+rc)o, and the singlet pairing amplitude, A r = 
>o, where r connects nearest (nn) or 



t J 



(cj t c 



i+rl 



il i+rf 



denotes the 



next nearest neighboring (nnn) sites; ( 
expectation value with respect to l^o)- 

This decoupling scheme of the renormalized Hamil- 
tonian (RMFT) leads to a BCS ground state |* ) = 

rifeK + «fc4t c -*i) with v t = i 1 - Ck/E k )/2 

and u| = 1 — v\. The corresponding gap equa- 
tions are £ r = — cos(k r)(k/ E/. and A r — 
l/Lj2k C0S (k r )Ak/Ek, together with the condition x = 
l/Lj^k (k/Ek for the hole-doping concentration 

Solving the gap equations, we find that a (i-wave pair- 
ing state is most stable for x < 0.4. In this case, 
A = |A X | = \A y \ with A x = -A y , £ = £ x = £„, and f = 
£x+y = £x—y The dispersion relation of the Gutzwiller- 
Bogoliubov quasiparticle is given by, Ek — y%| + Af , 
where, 



2g t i + J — xi + J' — X2 ) (cos fe^ + cos 
2<7 t t' + J" — xi + J' -x 2 ) 2 cos fen; cos k y 



A, = 



A 



(1) 



J — [3g s + 1 — (3 + x)ga] (cos k x — cos fc v ) . (2) 



In Eq. m and Eq. @ , the Gutzwiller factors correspond- 
ing to the kinetic, superexchange and the three-site term 
are g t = 2x/(l+x), g 8 = 4/(l + x) 2 and g 3 = Ax/(l + x) 2 
respectively Q. The last sum in Eq. Q is a sum 
over all pairs of neighboring sites r and r', where t T 
and t T i are the nn and nnn hopping terms. We define, 
J = 4t 2 /U, J' = it't/U, and J" = U /2 /U and abbre- 
viate, x\ = 3g s — 1 + 3(3 — x)gs, x 2 = 4(3 — x)ga, and 
X D = (1- x 2 )g 3 in Eq. (TJ. 

Fig. Ofa) shows the doping dependence of the mean 
field hopping parameters, £ and and the d-wave mean 
field pairing, A, obtained by solving the RMFT gap equa- 
tions. Since gt, 33, x 2 , xd and £' vanish as x — > 0, the 
dispersion relation, — > — (11/4) J£ (cos fe^ + cos fej/) in 
the limit of zero doping. 

The d-wave order parameter A decreases nearly lin- 
early with doping and vanishes around x 0.4, for our 
choice ofparameters. This is in agreement with previous 
studies 0[l3|- Note that we retain all terms of 0(t 2 /U) 
in the effective Hamiltonian, i? c fr, including the three-site 
term, which suppresses d-wave pairing in the overdoped 
regime. 

The doping dependence of the superconducting (SC) 



-4i4+^)\ Here, 
represents the expecta- 



gap, |Afe| at k = (jr, 0), is shown in Fig. {Rb). The 
doping dependence resembles experimental observations 
quite well. However, the gap is overestimated by a factor 
of about 2 [see scaling factor a = 1/2 in Fig-E^b)] within 
mean field theory which neglects additional off site cor- 
relations as well as dynamical effects due to the motion 
of holes 0| . The SC gap is not identical to the_true SC 
order parameter, $ = |{c||cj +T ^ 
t is a neighboring site and 
tion value calculated with the retransformed wave func- 
tion, exp(iS)P G \^o}- Calculating $ to 0{t/U), we find, 
$ m g t A + (t/U) g 3 (6 -x)A£, where we set t' « 0, for 
simplicity. As shown in Fig.QJb), ^ vanishes as x — > 0, 
as it should. 

We now consider the nature of the low lying excita- 
tions, the quasiparticles created at the nodal point, kp. 
The nodal dispersion around fop is characterized by the 
velocity, vp, which directly influences a number of ex- 
perimentally accessible quantities. Within RMFT, vf is 
directly obtained by calculating the gradient of (k along 
the direction, (0,0) — > (ir, n). The result is presented in 
Fig. 13a) (for t = 0.3, 0.4, 0.5 eV and a a = 4l) and is 
well approximated by the formula, 



Vp / o,q « v^sin/cF 



2g t (t + 2t' cos k F ) + XijZ 



(3) 



In the above equation, we set J', J", and 10 to zero for 
simplicity. As seen in Fig. !)p increases with x, but 
remains finite as x — ► 0. As can be inferred from Eq. J2J, 
the energy scale of the nodal velocity at x = is deter- 
mined by J, i.e., Vp/(aoJ) ~ \/2sin(fcp<)if£ w 1.5 (for 
£ w 0.38 and kp w ^). The observed doping dependence 
stems from the effects of Gutzwiller projection Pq. As 
x increases, holes gain kinetic energy by direct hopping, 
viz., g t increases with doping; but g s decreases, leading 
to the doping dependence of vp seen in Fig. [21 

Our results agree with the numerical VMC results of 
Paramekanti et ai, who extract vp from the discontinuity 
of the first moment of the spectral function in the repul- 
sive U Hubbard model |3j , and those of Yunoki et al. , who 




FIG. 1: (color online) (a) Doping dependence of the dimen- 
sionless mean field parameters £, A; (b) Doping depen- 
dence of (solid) the SC order parameter, $, and (dashed) the 
SC gap, a\A k \, at k = (tt,0) for t = 300 meV. The RMFT 
SC gap is scaled by a factor a = 1/2 for comparison with 
experimental data (red circles, Bi2212 Q]) 
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FIG. 2: (color online) Doping dependence of (a) Fermi ve- 
locity, vf, and (b) renormalization, Zy., of the Gutzwiller- 
Bogoliubov nodal quasiparticle. RMFT results are compared 
with experiments (in (a), data from Ref. Il9f) and VMC (in 
(b), data from Ref. 0), respectively. 



cited state in a projected superconductor is given by 



exp(i5)|4 r ^") = exp(iS)P'y_ k <\ i f?o), where the corre- 
sponding Bogoliubov QP operator is defined by -y_ fc , = 

WfccL fc ^ + UfcCfcf . In ARPES, the spectral weight corre- 
sponding to this excitation is determined by the matrix 
element MZ = \(^Z\c k ^)\ 2 /(NZN G ), where, NZ and 
Nq are the norms of \^Z) and I'F), respectively. Here, 
Cfe.cr = exp(— iS)ck,a exp(iS). Using the the Gutzwiller 
renormalization scheme we find, MZ ^ Z k n krT +O(t/U) 2 
[13 and, 



(4) 





for the QP renormalization [2l|, with e° = 2t(cosk x 



>(0,0) 



FIG. 3: (color online) (a) Energy dispersion, u — —E k , (b) 
quasiparticle weight, Zfen.^ CT , and (c) momentum distribution, 
{nka), of the Gutzwiller-Bogoliubov quasiparticle for different 
doping x; The energy dispersion [in (a)] is only shown when 
the corresponding quasiparticle weight is finite [see (b)]. The 
corresponding Fermi surface, ^ = 0, is shown in the inset of 
(a). 



obtain vp from the quasiparticle dispersion in the t — J 
model ,8] . A comparison to ARPES data |l9j , presented 
in Fig. |2f a), shows good agreement. The doping depen- 
dence of vf in the severely underdoped regime remains 
to be settled experimentally. While some groups report 
a nearly constant Fermi velocity (see data for LSCO in 
Fig. [2Jl, others observe an increase with doping (see data 
for YBCO and Bi2212 in Fig.0). Within RMFT, we also 
verify that the nodal properties remain essentially un- 
changed when A is set to 0; i.e., the doping dependence 
of vf results from the vicinity of the state to a Mott in- 
sulator, rather than the occurrence of superconductivity 
itself. 

We now calculate the spectral weight of the Gutzwiller- 
Bogoliubov quasiparticle (QP). The variational ex- 



vi, is the mo- 



cos k y ) + 4 1' cos k x cos k y . Here, n° CT 
mentum distribution in the unprojected wave function 
l^o) ■ The renormalization, Z k , of the nodal QP weight 
is plotted as a solid line in Fig.^b), and agrees well with 
VMC results The dashed lines correspond to results 
without i/L/-corrections. The dotted line, Z k — x, is the 
result from slave boson mean field theory (SBMFT). 

As a qualitative comparison with ARPES we show 
the energy dispersion, u> = —Ek, of the Gutzwiller- 
Bogoliubov QP along the directions, (0,0) — * (n,0), 
(ir, 0) — > (77, 7r), and (%, tt) —> (0,0) for different x in 
Fig. |3J We emphasize that our calculations describe the 
low energy sector, and do not seek to explain the "kink" 
at higher energies |l9j|. 

The spectral weight of the coherent peak, measured in 
ARPES, is related to the QP weight, M~ « Z k n° ka ; it is 
shown in Fig.|SJb). As seen in the figure, the QP spectral 
weight is severely modified by Gutzwiller projection. It 
decreases with doping, and vanishes at half filling. This 
causes a shift of spectral weight to an incoherent back- 
ground as seen in the momentum distribution function, 
(n ka ) fa Z k v\ + rif* + 0(t/U) 2 . While the first term 
corresponds to the coherent QP weight, the second gives 
the distribution of the incoherent part. We get, 



(1 - xf 
2(1+ a;) 

3 + x 



93- 



+ E 

T 



2U 



cos(£+ 
3g s - 1 



(1-x) 3 



35s 



33" 



(5) 



which is a smooth function of k. Results are shown in 
Fig-EI c ) ■ The incoherent weight is spread over the entire 
Brillouin zone, and overlies the coherent part from the 
Gutzwiller-Bogoliubov quasiparticles. At half-filling, all 
weight becomes incoherent. 

Finally, we consider the coherent QP weight, MZ ~ 
Z k v k , at the antinodal point, k = (7r,0). As seen in 
Fig. El» and Fig. |U it exhibits a non monotonic be- 
havior as a function of doping. Within our theory, this 
effect arises from a combination of the effects due to 
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FIG. 4: (color online) Doping dependence at the antinodal 
point, k = (7r,0): (a) QP renormalization, Zk, the unrenor- 
malized QP weight, n% a — v\, and the renormalized coherent 
QP weight, ZkV%; (b) coherent weight, Zuv\ compared with 
the experimentally determined Superconducting Peak Ratio 
(SPR) for Bi2212 Q. 



ior, when the underlying Fermi surface changes topology. 
Our results can be checked by experimental observations 
from ARPES of high temperature superconductors, thus 
providing a way to study the applicability of projected 
wave functions to these systems. The method we use is 
also amenable to other extensions such as coupling to the 
lattice, antiferromagnetism and long range interactions. 

We thank T. Valla, P. D. Johnson, J. Fink, and 
A. A. Kordyuk for discussions. VNM acknowledges par- 
tial support from the PSC-CUNY Research Award Pro- 
gram. BE thanks the Hermann Willkomm Stiftung for 
supporting his stay in Princeton. 



Gutzwiller projection and the topology change [see in- 
sert of Fig. Ola)] of the underlying Fermi surface (FS). 
Fig. Eta) illustrates this clearly. While the QP weight 
renormalization, Zk, increases with increasing doping, 
n° k = Vp., decreases due to the topology change, which 
occurs at x 0.15 — 0.20 for our choice of hopping pa- 
rameters (t' = —t/A). The change of the FS seems to be 
a generic feature of hole doped cuprates , although 
the exact doping concentration x, for which this occurs, 
is sensitive to the ratio between various hopping param- 
eters. The combined effect of strong correlations and 
topology change leads to a maximum of the QP weight for 
the doping level, x, at which the underlying FS changes 
topology. This result should be tested by appropriate 
analyses of recent ARPES data [22J . Indications for such 
a behavior have already been published |U 0] . Feng et 
al. [2^] extracted the Superconducting Peak Ratio [SPR, 
illustrated in Fig. EJb)] which is proportional to the co- 
herent QP spectral weight, Z^vf.. They found that the 
SPR increases with small x, attains a maximum value 
around x ~ 0.2 where it begins to decrease. Ding et al. 
[2^, reported similar results from ARPES. Although the 
topology change does not influence the stability of the SC 
state within RMFT, the SC pairing parameter $ (related 
to T c ) and the QP weight, Z^v\, show some similarity as 
a function of doping. 

To summarize, we studied the electronic structure of a 
Gutzwiller projected d-wave superconductor. A system- 
atic combination of renormalized mean field theory and 
direct evaluation of matrix elements was applied to the 
Hubbard model in the strong coupling limit. Our ana- 
lytical results can be used to fit experimentally observed 
quantities as well as those obtained from numerical meth- 
ods. The dispersion of the Gutzwiller projected super- 
conductor is renormalized in the vicinity of the Mott in- 
sulator, but the nodal Fermi velocity stays finite as the 
insulating limit is approached. The spectral weight of the 
nodal quasiparticle increases with doping, whereas that 
of the antinodal excitation shows non monotonic behav- 
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